An optical diode and magnifier from a general function photonic crystals 
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We have presented a new general function photonic crystals (GFPCs), which refractive index is a 
function of space position. Based on Fermat principle, we achieve the motion equations of light in 
one-dimensional general function photonic crystals, and calculate its transfer matrix. In this paper, 
we choose the line refractive index function for two mediums A and B, and obtain new results: (1) 
when the line function of refractive indexes is up or down, the transmissivity can be far larger or 
smaller than 1. (2) when the refractive indexes function increase or decrease at the direction of 
incident light, the light intensity should be magnified or weaken, which can be made light magnifier 
or attenuator. (3) The GFPCs can also be made optical diode. The new general function photonic 
crystals can be applied to design more optical instruments. 
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1. Introduction 

Photonic crystals are artificial materials with periodic 
variations in refractive index that are designed to affect 
the propagation of light [1-4]. An important feature of 
the photonic crystals is that there are allowed and for- 
bidden ranges of frequencies at which light propagates in 
the direction of index periodicity. Due to the forbidden 
frequency range, known as photonic band gap (PBG) [5- 
6], which forbids the radiation propagation in a specific 
range of frequencies. The existence of PBGs will lead to 
many interesting phenomena, e.g., modification of spon- 
taneous emission [7-9] and photon localization [10]. Thus 
numerous applications of photonic crystals have been 
proposed in improving the performance of optoelectronic 
and microwave devices such as high-efficiency semicon- 
ductor lasers, right emitting diodes, wave guides, optical 
filters, high-Q resonators, antennas, frequency-selective 
surface, optical limiters and amplifiers [11-18]. 

In Ref. [19], we have proposed special function pho- 
tonic crystals, which the medium refractive index is the 
function of space position, but the function value of re- 
fractive index is equal at two endpoints of every medium 
A and B. In this paper, we present a new general function 
photonic crystals (GFPCs), which refractive index is a 
arbitrary function of space position. Unlike conventional 
photonic crystals (PCs), which structure grow from two 
materials, A and B, with different dielectric constants ea 
and £b- Firstly, we give the motion equation of light in 
one-dimensional GFPCs based on Fermat principle. Sec- 
ondly, we calculate the transfer matrix for the GFPCs, 
which is different from the transfer matrix of the conven- 
tional PCs. Thirdly, we give the band gap structure and 
transmissivity. Finally, we choose the linearity refractive 
index functions for two medium A and B, and give the 
light field distribution in the GFPCs. We obtain some 
new results: (1) when the line function of refractive in- 



dexes is up, the transmissivity can be far larger than 1. 
(2) when the line function of refractive indexes is down, 
the transmissivity can be far smaller than 1. (3) when 
the refractive indexes function increase at the incident di- 
rection of fight, the light intensity should be magnified, 
which can be made light magnifier. (4) when the refrac- 
tive indexes function decrease at the incident direction 
of light, the light intensity should be weaken, which can 
be made light attenuator. (5) The GFPCs can be made 
optical diode, which transmits light from an input to an 
output, but not in reverse direction. 

2. The light motion equation in general func- 
tion photonic crystals 

For the general function photonic crystals, the medium 
refractive index is a periodic function of the space posi- 
tion, which can be written as n{z), n{x, z) and n(x, y, z) 
corresponding to one-dimensional, two-dimensional and 
three-dimensional function photonic crystals. In the fol- 
lowing, we shall deduce the light motion equations of the 
one-dimensional general function photonic crystals, i.e., 
the refractive index function is n — n{z), meanwhile mo- 
tion path is on xz plane. The incident light wave strikes 
plane interface point A, the curves AB and BC are the 
path of incident and reflected light respectively, and they 
are shown in FIG. 1. 

The light motion equation can be obtained by Fermat 
principle, it is 



S / n{z)ds = 0. 
J A 



(1) 



In the two-dimensional transmission space, the line ele- 
ment ds is 



ds = ^{dxY + {dzf = \J\ + Pdx, 



(2) 
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FIG. 1: The motion path of hght in the medium of refractive 
index n{z). 



where z = j^, then Eq. (1) becomes 

5 / n{z)^Jl + {zfdx = Q 

J A 

The Eq. (3) change into 



(3) 



dz 



dz 



At the two end points A and B, their variation is zero, 
i.e., Sz{A) = 5z{B) = 0. For arbitrary variation 5z, the 
Eq. (4) becomes 



dn{z) 
dz 

-n{z)- 



J2_^^2(^^ .2)- 

dz 

"I2_i2^(l + i2) 



simphfy Eq. (5), we have 

dn{z) zdz 



= 0, 



n{z) 



1 + z' 



(5) 



(6) 



The Eq. (6) is hght motion equation in one-dimensional 
function photonic crystals. 

3. The transfer matrix of one-dimensional gen- 
eral function photonic crystals 

In this section, we should calculate the transfer matrix 
of one-dimensional general function photonic crystals. In 
fact, there is the reflection and refraction of light at a 
plane siirfacc of two media with different dielectric prop- 
erties. The dynamic properties of the electric field and 
magnetic field are contained in the boundary conditions: 
normal components of D and B are continuous; tangen- 
tial components of E and H are continuous. We consider 
the electric field perpendicular to the plane of incidence, 
and the coordinate system and symbols as shown in FIG. 
2. 

On the two sides of interface I, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, there are 
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FIG. 2: The hght transmission and electric magnetic field 
distribution figure in FIG.l medium. 



Eq = Ej = Eti + E'^2 

Ho=Hi = Hn cos 6^ - H'^^ cos 61 ^ ' 



On the two sides of interface II, the tangential compo- 
nents of electric field E and magnetic field H are contin- 
uous, and give 



Ell = E'j = Ei2 + Er2 



(8) 



Hii = H'j = Hi2 cos 6»/ - Hr2 COS 0/ , 

the electric field En is 

Ea = ^tioe*(*=»^^+'=^^)U=o = ^tioe' ^"(o)-"^^-^, (9) 
and the electric field Ei2 is 

^ ^^^i^n{b){sinOlxB-\-cos0fb) (10) 

Where xa and xb are x component coordinates corre- 
sponding to point A and point B. We should give the 
relation between Ei2 and En- By integrating the two 
sides of Eq. (6), we can obtain the coordinate compo- 
nent Xb of point B 

/■"(^^ dn{z) _ r''' zdz 
Jn{0) n{z) l-hi^' 



to get 



and 



dx 



dz 



(12) 



(13) 
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FIG. 3: The structure {BA)^ of the general function photoni( 
crystals. 



where fco = cot^^ and A:^ = ^ From Eq. (12), there i; 
n{z) > n(0) sin $1 . and the coordinate ajg is 



Xb=Xa + 



f 

Jo 



dz 



(14 
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where b is the medium thickness of FIG. 1 and FIG. 2. 
By substituting Eqs. (9) and (14)into (10), and using th( 
equality 



we have 



where 



n{0)smej. = n{b) sin6lj, 



Ei2 = Ene'^\ 



Sb = -nb{b) (cos 0l b + sin 9^ i 
c Jo 

and similarly 



dz 



- 1 



E' o = Er2e 



(15j 
(16) 

= ), (17) 
(18) 



Substituting Eqs. (16) and (18) into (7) and (8), and 



using H = ./j^nE, we obtain 



where 
Mb = 



cosdh 



(19) 



^nb(b) cos 0' ,, , 
^0 ' • 1(^0) 



-m,(0)yi^cos^,^sin5, 



The Eq. (20) is the transfer matrix M in the medium of 
FIG. 1 and FIG. 2. By refraction law, we can obtain 



smf 



no 

— — smc- jcost/j = 4/ 1 - - „ 

«,(0) ' V (0) 



sin 61°, cosel = \ 1- sin^ 0^ (21) 



where no is air refractive index, and n(0) = n{z)\z=o- 
Using Eqs. (15) and (21), we can calculate cos9f. 

4. The transmissivity cind light field distribu- 
tion of one-dimensional general function photonic 
crystals 
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FIG. 4: The two periods transmission figure of light in general 
function photonic crystals. 



In section 3, wo obtain the M matrix of the half pe- 
riod. We know that the conventional photonic crystals 
is constituted by two different refractive index medium, 
and the refractive indexes arc not continuous on the in- 
terface of the two mediums. We could devise the one- 
dimensional general function photonic crystals structure 
as follows: in the first half period, the refractive index 
distributing function of medium B is nb{z). and in the 
second half period, the refractive index distributing func- 
tion of medium A is na(z), corresponding thicknesses are 
b and a, respectively. Their refractive indexes satisfy con- 
dition nb{h) 7^ na(0), their structure are shown in FIG. 
3, and FIG. 4. The Eq. (20) is the half period transfer 
matrbc of medium B. Obviously, the half period transfer 
matrix of medium A is 
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FIG. 5: The line refractive index functions in a period. The 
FIG .5(a) is the up line function of refractive indexes, and 
FIG .5(b) is the down line function of refractive indexes. 
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In one period, the transfer matrix M is 
M = Mb • Ma 



^ cos 6( sin Sh 




rib 



— i sin (5b 
/^nt{b) cos 
,(0) cos 0^ cos Sj, 
nb{b) cos Of 

— i sin Sg 
/ ^na(a) cos U-- 

a (0) COS COS (5t 

7la (a) COS 9^^ 




The form of the GFPCs transfer matrix M is more com- 
plex than the conventional PCs. The angle O^, Oj, 
and are shown in Fig. 4. The characteristic equation 

of GFPCs is 
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Hi 
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En+1 
Hn+1 
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En+1 
Hn+1 



(28) 



Where N is the period number. With the transfer 
matrix M (Eq. (28)), we can obtain the transmission 
and reflection cocfFiciGnt t and t, and the transmissivity 
and reflectivity T and R, they are 
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(29) 
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FIG. 6: The relation between transmissivity and frequency 
corresponding to the up line function of refractive indexes 
(FIG .5(a)) 



and 



T = t-t*, 
R = r-r*. 



(31) 
(32) 

In the following, we 

give the electric field distribution of light in the one- 
dimensional GFPCs. The propagation figure of light in 
one-dimensional GFPCs is shown in FIG. 9. From Eq. 
(28), we have 



Where rjo = t]n+i = 



El 
Hi 



= Mi{di)M2{d2) ■ ■ ■ MN-iidN-i) 
EN{di 



1(33) 



/A \l J^N{U,i-td2 \-dN-l+AzN) 

Mn{Azn) jj^^^^ +d2---+ dN-l + Azn) 

where di and ^2 arc the thickness of first and second 
period, respectively, Azn is the propagation distance of 
light in the N-th period, Ei and Hi are the intensity of 
incident cilcictric field and magnetic field, and EN{di + 

d2 h djv-i + A^;;v) and HN{di +d2 h <ijv-i + Azn) 

are the intensity of the N-th period electric field and mag- 
netic field. The Eq. (34) can be written as 



EN{di + d2- ■ ■ + dN-i + Azn) 
HN{di +d2 1- djv-i + Azn) 



M^\Azn) 



M^\{dN-i) ■ ■ ■ M^\d2)M^\di) { 



( A{Azn) B{Azn) 
\C{Azn) D{Azn) 



El 
Hi 



(34) 



the electric field Ei and magnetic field Hi can be written 
as 



Ei=Eii+Eri = {l + r)Eii, 



(35) 



5 



Hi = Hii cos 6° - Hri cos 

= .[^cose1{l-r)Ea. (36) 
V Mo 

From Eqs. (34)-(36), we can obtain the ratio of the elec- 
tric field EN{di+d2- • •-|-djv-i + A^jv) within the GFPCs 
to the incident electric field En, it is 

I E^idi + C?2 h + ^Zn) |2 

En 

= \A{Azn){1 +r) + B{AzN).f^cos9°{l - r)|^7) 

V Mo 

5. Numerical result 

In this section, we report our numerical results of trans- 
missivity. We consider refractive indexes of the linearity 
functions in a period, it is 

n(,(^) = n6(0) + ^^^fci^^, 0<z<b, (38) 

n„(^) = na(0) + ^°^^^~""^°^ z, 0<z<a, (39) 
a 

Eqs. (38) and (39) are the line refractive indexes distri- 
bution functions of two half period mediums B and A. 
When the Midpoint values nb(0), ribih), na(0) and na{a) 
are all given, the line refractive index functions n},{z) and 
na{z) are ascertained. The main parameters are: the half 
period thickness b and a, the starting point refractive in- 
dexes nb(0) and na(0), and end point refractive indexes 
rH){h) and na{a), the optical thickness of the two medi- 
ums are equal, i.e., 715(0)6 = na(0)a, the incident angle 
6^ = 0, the center frequency ujq = 1.215 x IQ^^Hz, the 
thickness b = 280nm, a = 165nm and the period number 
= 16. 

In FIG. 5(a), we take nb(0) = 1.38, rn,{b) = 1.9 for 
the medium B, and na(0) = 2.35, na{a) = 4.25 for the 
medium A, which are the up line function of refractive 
indexes. In FIG. 5(b), we take nh{Q) = 4.25, nb{b) = 2.35 
for the medium _B, and na(0) = 1.9, na{a) ~ 1.38 for the 
medium A, which are the down line function of refractive 
indexes. By the refractive indexes function, we can cal- 
culate the transmissivity. With the FIG. 5(a) and 5(b), 
we obtain the transmissivity distribution in FIG. (6) and 
FIG. (7). From FIG. (6) and FIG. (7), we obtain the re- 
sults: (1) when the line function of refractive indexes is 
up, the transmissivity can be far larger than 1 (T max- 
imum is 10^). (2) when the line function of refractive 
indexes is down, the transmissivity can be far smaller 
than 1 (T maximum is lO"'^). In the following, we shall 
study the light field distribution of the one-dimensional 
GFPCs for the light of positive and negative incident. 
The positive incident is shown in FIG. 9 and the nega- 
tive incident is shown in FIG. 10. The refractive indexes 
line function of positive incident is in FIG. 5(a), and then 
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FIG. 7: The relation between transmissivity and frequency 
corresponding to the down line function of refractive indexes 
(FIG .5(b)). 
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FIG. 8: The input and output light in the GFPCs. 



the refractive indexes line function of negative incident 
is in FIG. 5(b). The FIG. 5(a) is the function of fine 
increasing, and the FIG. 5(b) is the function of line de- 
creasing. From Eq. (37), we can calculate the electric 
field distribution of light in the GFPCs. For the posi- 
tive incident, the electric field distribution is shown in 
FIG. 11(a), and FIG. 11(b) is the electric field distribu- 
tion for the negative incident. From FIG. 11(a) and (b), 
we can obtain the result: (1) when the refractive indexes 
function increase at the incident direction of light, the 
light intensity should be boosted up or magnified, which 
can be made light magnifier (magnifying multiple 10^). 
(2) when the refractive indexes function decrease at the 
incident direction of light, the light intensity should be 
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FIG. 9: The light positive incident to the GFPCs. 
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FIG. 10: The light negative incident to the GFPCs. 
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FIG. 11: The light distribution in the GFPCs. Figure (a) is 

corresponding to the positive incident (FIG. 9), and Figure 
(b) is corresponding to the positive incident (FIG. 10). 



■weaken, 'which can be made light attenuator (magnifying 
multiple less than I). Actually, the GFPCs structure in 
FIG. 10 is the optical diode, since the light intensity be 
magnified as the light positive incident, and the light in- 
tensity should be decreased as the light negative incident. 

6. Conclusion 

In summary, We have theoretically investigated a 
ne-w general function photonic crystals (GFPCs), 'which 
refractive index is a function of space position. Based 
on Fermat principle, 'we achieve the motion equations 
of light in one-dimensional general function photonic 
crystals, and calculate its transfer matrix. We choose 
the line refractive index function for two mediums A and 
B, and obtain some results: (1) when the line function 
of refractive indexes is up, the transmissivity can be far 
larger than 1. (2) 'when the line function of refractive 
indexes is do'wn, the transmissivity can be far smaller 
than 1. (3) when the refractive indexes function increase 
at the direction of incident light, the light intensity 
should be magnified, 'which can be made light magnifier. 
(4) when the refractive indexes function decrease at the 
direction of incident light, the light intensity should 
be weaken, which can be made light attenuator. (5) 
The GFPCs can be made light diode. The new general 
function photonic crystals can be applied to design more 
optical instruments. 
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